Mott-insulator phase of coupled ID atomic gases in a 2D optical lattice 
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We discuss the 2D Mott insulator (MI) state of a 2D array of coupled finite size ID Bose gases. 
It is shown that the momentum distribution in the lattice plane is very sensitive to the interaction 
regime in the ID tubes. In particular, we find that the disappearance of the interference pattern in 
time of flight experiments will not be a signature of the MI phase, but a clear consequence of the 
strongly interacting Tonks-Girardeau regime along the tubes. 
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Remarkable developments in atomic cooling and trap- 
ping have triggered the interest in strongly correlated 
atomic systems 0,0,0,0,0- In particular, the physics of 
cold atoms in periodic potentials induced by laser stand- 
ing waves (optical lattices) has attracted a major atten- 
tion, mostly due to its links to solid state physics , and 
due to the observation of the superfluid (SF) to Mott- 
insulator (MI) transition Q in Munich experiments ||. 

The reduction of spatial dimensionality in these sys- 
tems is now a "hot topic", in particular with regard to 
the creation of ID gases, where the interaction between 
particles becomes more important with decreasing the 
gas density. For low densities or large repulsive interac- 
tions the system enters the strongly-interacting Tonks- 
Girardeau (TG) regime, in which the bosons acquire 
fermionic characteristics and dynamical and correla- 
tion properties drastically change 
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regime requires tight transverse trapping, low atom num- 
bers, and possibly the enhancement of interactions via 
Feshbach resonances 0, El In this sense, 2D op- 

tical lattices are favorable, since the on-site transverse 
confinement can be made very strong and for sufficiently 
small tunneling rate each lattice site behaves as an in- 
dependent ID system. Recent experiments on strongly 
correlated ID g ases have been performed along these lines 

These studies motivate the analysis of an interesting 
physics in a (2D) array of coupled ID Bose gases. In 
a 2D lattice the coupling is provided by the inter-site 
tunneling, and each lattice site is a ID tube filled with 
bosonic atoms. This regime is easily achievable experi- 
mentally by lowering the lattice potential, and it repre- 
sents the bosonic analog of ID coupled nanostructures 
[lflj . As was first shown by Efetov and Larkin ji^J], for 
infinitely long ID tubes at zero temperature any infinites- 
imally small tunneling drives the system into the super- 
fluid phase. The gas then enters an interesting cross- 
dimensional regime, in which it presents ID properties in 
a 3D environment |2 lL |22| . For ID tubes of finite length 
L, at sufficiently small tunneling rate the system can un- 
dergo a cross-over from such anisotropic 3D superfluid 



state to the 2D Mott insulator state |22j . Strictly speak- 
ing, this 2D MI phase requires a commensurable filling of 
the tubes, i.e. an integer average number of particles N 
per tube. Then the system of finite-length tubes at zero 
temperature is analogous to that of infinite tubes at a 
finite temperature T, and the critical tunneling t c for the 
T = cross-over to the MI phase can be obtained from 
the finite-temperature results of Ref. [2(J by making a 
substitution 1/T -> L [H- 

This Letter is dedicated to the analysis of correlation 
properties of this 2D Mott insulator. We show that the 
momentum distribution is crucially modified by a com- 
bined effect of correlations along the ID tubes and inter- 
tube hopping. For the case of a weakly interacting gas in 
the tubes, the phase coherence is maintained well inside 
the MI phase. This is similar to the situation in 2D and 
3D lattices, studied by means of Quantum Monte Carlo 
calculations |2^j and investigated experimentally through 
the observation of an interference pattern after switching 
off the confining potential |||. However, an increase of 
the interaction between particles in ID tubes reduces the 
inter-tube phase coherence and flattens the momentum 
distribution in the transverse direction(s). In particular, 
the interference pattern observed in Ref. || should be 
largely blurred if the ID tubes are in the TG regime. 
This effect can be revealed in current experiments. Re- 
markably, the disappearance of the interference pattern 
will not be a signature of the MI phase, but a clear con- 
sequence of the strongly interacting regime for ID tubes. 

In the following we consider a Bose gas at zero temper- 
ature in a 2D optical lattice, such that every lattice site 
can be considered as an axially homogeneous ID tube of 
finite size L, with N = nL being the number of parti- 
cles per tube, and n the ID density. The tunneling be- 
tween neighboring tubes is characterized by the hopping 
t, which depends on a particular lattice potential and 
atomic species employed. We label tubes by the index j 
and denote by x the coordinate along the tubes. 

The action describing the coupled tubes has the form: 
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/oo fL/2 
dr dx (^j+^i), (1) 
-oo 7-i/2 

where ipj(x,T),ij)j(x,T) are complex bosonic fields asso- 
ciated with the j-th tube, the symbol < ij > denotes 
nearest neighbors, and Sj is the action describing the 
physics along the j-th tube. 

In the absence of tunneling, there are no correlations 
between different tubes, and the one-body Green function 
is diagonal: 

Gij{xi-X2,T1-T%) = (lpi{xi,Ti)lpj(x2,T2)} 

= 6ijGo(xi-X2,Tl-T2). (2) 

The presence of tunneling between neighboring tubes, 
provided by the second term on the rhs of Eq. l[T]l. modi- 
fies the momentum distribution. Above a critical tunnel- 
ing amplitude t c the system undergoes a cross-over from 
the MI to an anisotropic 3D SF phase j^- This cross- 
over and the MI phase can be analyzed within the ran- 
dom phase approximation (RPA) 1251. successfully used 
in the studies of coupled spin chains [2g] . Decoupling the 
tunneling term in the action © by using the Hubbard- 
Stratonovich transformation and keeping only the leading 
quadratic terms, yields the RPA Green function in the 
momentum- frequency representation: 



G(q,k,co) 



G (k,w) 



l-T(q)G (k,u) 



(3) 



with q = (q y , q z ) being the quasimomentum in the lattice 
plane, T(q) = 2t(cosq y a + cosq z a) 1 a the lattice constant, 
and Go(k, u>) the Fourier transform of the Green function 
(J2J) (hereinafter we put h = 1): 



G (fc,w) 



r L/2 

dr I dx e 

-L/2 



G (x,t). (4) 



The long-wavelength behavior of the Green function 
Gq{x, t) can be found using Luttinger liquid theory [27| . 
At zero temperature, employing a conformal transforma- 
tion in order to take into account the finite size L of the 
tubes 25], we obtain: 



G (x,t) = n 



n 2 /N 2 



sinh (tt(/L) sinh (ir(/L) 



(5) 



where £ = v s t + ix, and v s is the sound velocity. The 
interactions enter Eq.JSJ) through the factor d — 1/AK 
related to the interaction-dependent Luttinger parameter 
K. The Fourier transform of Eq. © yields 



1 / N 
G (fc,c) = — - 
nv s \ Zir 



2-2d 



ujL 



2ir 2irv* 



(6) 



where the quantity I(p, fl) is expressed through the hy- 
pergeometric function 3^: 



Re 



P i r(d) 

/ [d.d-p.'^r^ 



d + p — iQ 



, (7) 



with p = kL/2ir and Q = luL/2ttv s being the dimension- 
less momentum and frequency. Integrating Eq. 10 over 
u) one obtains the momentum distribution No(kL /2tt) = 
J du>Go{k, lu)/2tt in the absence of tunneling: 



N 1 2tt / 



-2d 



pW(d) 



2 F 1 (d,d + p;l+p;l), (8) 



which behaves as p 2 ^ 1 for p > 1. The Luttinger liq- 
uid description employed here is valid for low momenta 
k -C mi. Accordingly, the dimensionless axial momen- 
tum p = fcL/27r, which is an integer number, should sat- 
isfy the inequality p <C N. The momentum distribution 
|JHJ represents the fraction of particles in the state with 
momentum p and is normalized as ^ p No(p) = N. 

The critical tunneling t c for the MI to SF cross-over 
is obtained as the value of t for which the denominator 
of Eq. @ vanishes for zero momenta k and q and zero 
frequency u). We thus have 



tr 



fi 4/i \2tt 



N 



2d-2 



1 



7(0,0) 



(9) 



Note that with t c from Eq. ©, the Green function G 
in Eq. (|SJ) becomes a universal function of the dimen- 
sionless quantities t/t c , p, qa and fi. For the TG regime 
of ID bosons in the tubes, the Luttinger parameter is 
K = 1 and d = 1/4. Then, as the chemical potential 
is fi = mv 2 = ir 2 n 2 /2m, from Eq.© we obtain i c //i ~ 
0.05A~ 3 / 2 . For the weakly interacting regime, the Lut- 
tinger parameter in the ID tubes is K = n^n/mg) 1 ^ 2 1 
and Eq.(0) gives 1(0, 0) = 16nK > 1. In this regime the 



chemical potential is /i 



ng, and Eq.© then 



yields i c //i ~ (1/16) A -2 (as expected from the mean- 
field calculations for a 2D lattice of zero-dimensional sites 
These results are in qualitative agreement with the 
recent calculations of Ho et al. [13 • One clearly sees that 
strong correlations along the tubes drastically shift the 
boundaries of the MI phase. 

The momentum distribution for the coupled ID tubes 
in the MI phase, N(q, k), is obtained by integrating the 
Green function © over the frequency. Our calculations 
show that only the lowest axial mode for which the mo- 
mentum k = 0, is significantly affected by the tunneling. 
The physical reason is that the k = mode is approach- 
ing the instability on approach to the critical tunneling t c , 
whereas k ^ modes are still far from instability. This is 
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Figure 1: Transverse momentum distribution for k = at 
t/t c = 0.3 for K = 1 (solid), K = 4 (dotted), and if = 25 
(dashed). In the figure we have chosen q y = q z = q. 



0.9 



H 0.8 



0.7, 



1 1 1 1 1 1 1 




s '■• ** ^ 


(a) 


N '*. * _ 

N 

X 

N 

N 

N 

S. 

N 



0.25 



0.5 



0.75 




Figure 2: Transverse momentum distribution for t/t c = 0.3, 
at K = 1 (solid), K = 4 (dotted), and K = 25 (dashed), for 
TV = 50 (a) and iV = 500 (b). In the figures we have chosen 

q y = q* = q- 



reflected in the resonance character of the Green function 

for fc = and t — > f c . Our results in Fig.Q]show that 
the transverse distribution N(q, 0) corresponding to the 
fc = axial mode is non-flat for any interaction regime 
along the tubes. In contrast, for k ^ the quantity 
T(q)Go(k,u) is always small. Therefore, expanding the 
rhs of Eq. @ in powers of T(q) up to linear order and in- 
tegrating over to we obtain an almost flat transverse mo- 
mentum distribution for k ^ modes: N(q,p)/N(0,p) = 

1 - {Ad 2 /4p 2 - d )(t/t c )(2 - Y,i= y , z coaq.a), where the co- 
efficient A is of order unity and the second (q-dependent) 



term is always very small. 

We now turn to the discussion of the transverse quasi- 
momentum distribution N±(q) = J2 k N(q, k). The sum- 
mation over the axial modes changes the picture dras- 
tically compared to the distribution for a given k. As 
only the k = component is significantly affected by the 
tunneling, one can rewrite Eq. © in the form: 



G{q,k,u)~G {k,u) 



l-T(q)G o (0,u 



Sk,o- (10) 



In the second term on the rhs of Eq. l|lUf) we may use 
the Green function Go(Q,uj) following from Eqs. (JSJ and 
(J7J). For fc = (p = 0), one can put the hypergeomet- 
ric function 3F2 = 1 in Eq.Q, which gives 1(0, £1) ~ 
4ird/(cP + O 2 ). Omitted terms give a very small rela- 
tive correction of the order of d 3 < 1/4 3 . Hence, us- 
ing Eq.@, for the Green function at fc = in the ab- 
sence of tunneling we have Gq(0,u>) — d 2 /4t c (d 2 + fi 2 ). 
Then, integrating Ea. i|l(J|) over f2, summing over the ax- 
ial modes fc, and imposing the normalization condition 
N = J^k I dujGo(k,oj)/2ir for the first term on the rhs, 
we obtain the transverse momentum distribution 



N±(q) 
N 



= 1+ ^ 



2d 




(11) 



normalized by the condition (a/27t) 2 / d 2 qN±(q) = N. 

In Fig. |2 we depict the results of Eq. l|llfl for different 
values of ./V and the Luttinger parameter K . Due to the 
prefactor in the second term on the rhs of Eq. the 
transverse momentum distribution strongly depends on 
the interaction regime along the tubes. 

For the weakly interacting regime (d <C 1), the dis- 
tribution N±(q) is not flat even deeply inside the MI 
phase. Similar results have been obtained by means of 
Quantum Monte Carlo calculations [531 f° r the case of 
lattices of zero-dimensional sites. In our case, only for 
rather low tunneling (t/t c < 0.1) the quasimomentum 
distribution becomes flat, and switching off the lattice 
potential should lead to a blurred picture as that ob- 
served by Greiner et al. 0. Non-flat distributions in 
the MI phase as those of Fig. [3 will manifest themselves 
through the appearance of interference peaks in the same 
type of experiment. 

On approach to the TG regime, the quasimomentum 
distribution becomes progressively flatter. The main rea- 
son for this behavior is that when the system becomes 
more interacting, the fc = component is more depleted, 
contributing less to the total quasimomentum distribu- 
tion. Therefore, if the ID tubes approach the strongly 
interacting TG regime, in experiments as those of Ref. [8| 
the interference pattern will be essentially smeared out. 
One may expect a partial destruction of the interference 
pattern even for moderate values of the Luttinger param- 
eter (see, e.g., the case K = 4 in Fig.|2J|. 
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The random phase approximation used in our calcula- 
tions, was shown to be a good approximation for a wide 
range of parameters of coupled one-dimensional Heisen- 
berg spin chains j2(| . Here we give yet another estimate 
for the applicability of RPA, relying on the Ginzburg 
criterion adapted to a quantum phase transition at zero 
temperature. We compare fluctuations of the order pa- 
rameter in a volume determined by the correlation radius 
extracted from the Green function 0, with the scale on 
which the non- linear effects become important. The lat- 
ter is obtained from the four-point correlation function 
of each tube. We have found that RPA is adequate for 
(t c —t)/t c 3> B(d), where B(d) has been obtained numeri- 
cally from the four-point correlation function. In the case 
of the Tonks-Girardeau regime along the tubes, we have 
B « 0.1, and it decreases significantly with decreasing d 
and entering the Gross-Pitaevskii regime. 

In conclusion, we have considered a bosonic gas in a 
2D optical lattice of finite ID tubes at zero temperature, 
focusing our attention on the momentum distribution in 
the MI phase. We have shown that the strong correla- 
tions along the ID tubes significantly modify the quasi- 
momentum distribution in the lattice plane. We have 
found that in the MI regime only the lowest momentum 
along the tubes is affected by the inter-site hopping, and 
hence only this component contributes to the formation 
of interference fringes. Consequently, the larger the in- 
teractions are (larger depletion) the less pronounced is 
the visibility of the interference fringes. In particular, 
for the TG regime in the tubes, the quasimomentum dis- 
tribution becomes progressively flatter, leading to an ob- 
servable blurring of the interference pattern after expan- 
sion. This effect can be observed in current time of flight 
experiments, and can be used to reveal a clear signature 
of the strong correlations along the sites. 
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